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Abstract 

We present the Generalized Borel Transform (GBT). This new approach al- 
lows one to obtain approximate solutions of Laplace/Mellin transform valid in 
both, perturbative and non perturbative regimes. We compare the results pro- 
vided by the GBT for a solvable model of quantum mechanics with those provided 
by standard techniques, as the conventional Borel sum, or its modified versions. 
We found that our approach is very efficient for obtaining both the low and the 
high energy behavior of the model. 
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Recently, Penin and Pivovarov M presented a numerical analysis of renormalon 
techniques in quantum mechanics. They used a simple solvable model, namely a delta 
function scattering potential, to address the problem of resummation of perturbative 
series @. In fact, this quantum mechanics potential can be considered as a confining 
one, mimicking the long sought property of Quantum Chromodynamics (QCD). This 
is one of the many examples where the finite order perturbation theory predictions 
present uncertainties comparable with experimental errors |§. In all these cases one 
is forced to take into account, in a way or another, non-perturbative contributions. In 
connection with these difficulties, it is certainly instructive to study exactly solvable 
models so that the efficiency and precision of the different proposals Jl|[|] can be 
quantitatively checked. This analysis can help us to define criteria for selecting the 
appropriate approach to be used in more realistic QCD is. 

In reference |l| the delta function model is tackled by means of standard Borel 
summation techniques combined with the renormalon approach |J. On the other hand, 
we have recently introduced || a Generalized Borel Transform (GBT) that avoids the 
implementation of a perturbative expansion. This proposal was successfully applied |F] 
to obtain an analytic expression for the heavy quark- ant iquark potential, valid for all 
distances. 

The purpose of this paper is to discuss the use and advantages of the GBT in the 
case of quantum mechanics and in particular in connection with the model studied 
in ref. [|TJ. We are able to obtain analytic expressions valid in all the range of the 
variables. 

We start by briefly summarizing ref. U in order to prepare the basic material of 
the delta potential that is further treated with our GBT approach. 

The quantum mechanics model, analyzed first in the scattering region, is defined 
through the singular potential 

V (r) = v5 (r — a) ; v, a = cte (1) 

and can be considered as a kind of confining (not completely) interaction. Usually, the 
dispersion relation related to this kind of singular potential is analyzed in analogy with 
dispersion relations in elementary particle theory Here, to simplify the presentation, 
we consider only the s - wave amplitude and the study of the wave function at the origin. 
To analyze the scattering of a wave packet, one considers an integral of the form 

/•OO 

*(Q) = 1 + F(Q)= / ^(q)W(q,Q)dq (2) 

Jo 

where, 

q 2 

is a momentum distribution function normalized to one and, 
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1 + - 

2q 



exp (— 2qa)) 



-i 



(3) 



is the exact solution for the scattering of a plane wave in the " Euclidean" region where 
the momentum q > is obtained through a Wick rotation. This is possible whenever 
m = l/a > \v | and consequently there are no bound states present. 

The perturbative solution in eq. (Q) suffers of the same difficulties present in per- 
turbative QCD and it mimics some general features of renormalons. 

We start the analysis with the perturbative treatment in order to make explicit the 
difficulties of the standard method. At high energies (q ^> m) one can write 



(4) 



n=0 



where a (q) = v/2q is the natural parameter of expansion in expression (|3]). 
Replacing now the series (|]) into (fj) one obtains 



tf« (g) = X>! (-«(?))' 



(5) 



n=0 



For v > , the theory is Borel summable || and the conventional Borel transform is 



°° 1 1 POO 

b (s) = (-*) n = — t - n s (?) = -m / 5 (s) exp {q)) (6) 

^ s + 1 a(g)io 

Consequently, the approximate solution (^) results 

*ST (?) = exp (l/a (g)) Ei (1, l/a (q)) (7) 

where [Q] 

Ei(^)^ eXp( ; to)dt ,= 1,2.. *>0 (8) 

A numerical comparison of the approximate expression (0) with the exact solu- 
tion presented in Figure 1, clearly shows that it is not a good approximation for 
intermediate values of q ~ m. 

An alternative approach |K| uses a modified perturbation theory and provides 
results of better precision at very small momenta. This technique starts by choosing, 
instead of a(q), a modified expansion parameter of the form 

a » (?) = 77— 1 



2q + fi \ 2q + fi 

where \l — fi (m, v) is fitted using experimental information. In this model, the exact 
solution ^ (Q) given in eq. @, plays the role of experiment data. 
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The accuracy of the approximation and the optimal value of [l clearly depend on 
the range where data have to be fitted. Unfortunately, the value that optimize the 
result at very small Q, does not fulfil this requirement at very large Q as it is evident 
in Figure 1. 

Hence, the method and its modification work reasonably well in providing an ap- 
proximation to the exact solution but each one in different regimes of Q. They are not 
able to reproduce simultaneously, with sufficient precision, the exact result in both IR 
and UV regimes. 

We have found that the Generalized Borel Transform f|[[7| copes with that diffi- 
culty. This GBT approach was designed to obtain the approximate solution of any 
Laplace/Mellin transform in the wide range of the parameter involved. 

In fact, when using our proposal, one is performing a whole class of transformations, 
depending on a parameter that we call A. As it should be, the result does not explicitly 
depends on A. Consequently, one can choose for this parameter the best adapted value 
in each particular problem under consideration. In practice, as one is performing an 
approximate calculation, one can get rid of A by means of a suitable saddle-point like 
technique. 

Let us briefly present our technology in connection with a Laplace-Mellin transform 
presents when treating the quantum mechanical problem under analysis. In this case, 
we typically face a transformation like 

POO 

S (g, a,n) = / x n H [x, a) exp (—gx) dx ; g > (9) 
Jo 

where we have explicitly extracted a factor x n from the function to be transformed, 
because, at it will became clear after, this operation facilitates the implementation of 
our proposal. 

The Generalized Borel Transform (GBT) of 5* is defined as 



B x (s,a,n) = J exp [s/rj\ 



~1 -As 



Xrj 



+ 1 



S{g,a,n)d{l/rj) ; Re (s) < (10) 



A being any real positive non zero value and where we have defined 

1/77 = A [exp (g/X) - 1] 



(11) 



This particular election of 77 allows one to define the function u\ (g) monotonically 
increasing in the interval < u\ (g) < 00, namely 



(g) = - ~ A In 
V 



Xrj ~^ 



(12) 



Consequently B\(s) results in 
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B\ (s, a, n) = y exp (su) S [g\ (u) , a, n] [1 + A r\ (u)] du ; Re (s) < 

Q 

The last expression can be written as a Laplace transform 

oo 

B\ (s, a,n) — J exp (su) L x [g (u) , a, n] du ; Re (s) < 
o 

of the function Lx [g (u) , a, n] there implicitly defined. 

Replacing (|9|) and ( pTTj) in equation (|T0|) we can easily test that the GBT is an 
analytic function on the negative Borel half-plane, such that its extension to the other 
half-plane also exists and is also analytic with a cut on the real positive axis. From 
this observation, S (g,a, n) can be expressed in terms of the inverse Laplace transform 
integrated on the above mentioned cut 

oo 

S (g, a, n) = ^ - J exp [su\ (g)] AB X (s, a, n) ds (13) 

o 

The main advantages of this proposal comes from the fact that it allows one to 
perform the calculations in terms of the parameter A. Moreover, as it will become 
clear below, this approach avoids the implementation of a perturbative expansion. Each 
value of the parameter A defines a particular Borel Transform. This can be summarize 
by writing S (g, a, n) = T^ 1 [T A (S (g, a, n))] where T A (S (g, a, n)) = B x (s, a, n) . 

The discontinuity of the Bx (s, a, n) can be expressed as 



A£> A (s, a, n) = 2A J dwexp[T (w,X,s,a,n) 



where 



T (w, A, s, a,n) = — In {T [A (s + x (w))]} + {A [s + x (w)] — 1} In (As) — As 
+w + In [x n H (x, a)] 

with x = exp (w) and F represents de Euler gamma function. After the change of 
variables s = A exp (t) one gets 

oo oo 

S (g, a, n) = 2A 2 (1 — exp (— g/X)) J J exp [G (w,t, g, A, a,n)} dwdt (14) 

— oo — oo 

with 
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G (w, t, g, A, a, n) = —s (t) u\ (g) + 1 + T [w, A, s (t) , a, n] 



(15) 



In this way we can obtain the dominant contribution of the double integral using 
the steepest descent technique in the combined variables [t,w]. In so doing, one first 
computes the saddle point t (g, a, n) and w Q (g, a, n) in the limit A ^> 1. In this case, 
the saddle point is 



t Q (g, a, n) = In 



x 2 (g,a,n) 



f [x (g,a,n),a,ri[ 



w (g,a,n) = In [x {g,a,n)] 



(16) 



where x Q (g, a, n) is the real and positive solution of the implicit equation coming from 
the extremes of the function G in the asymptotic limit in A, namely 



2 2 

x Q g 



and 



f(x ,a,n) [f {x ,a,n) + 1] 
din [H (x , a)} 



f (x 0} a } n) = 1 + n + x 



(17) 



(18) 



In the range of the parameters where / [x Q (g)] 3> 1 , we have retained the first 
order in the expansion of G around the the saddle point. In so doing, it remains 



S Ap (g,a,n)^4n\ 2 [1 - exp {-g/X)\ 



exp [G (w ,t ,g,\,a,n)] 
\J D (t a , w ) 



(19) 



where 



D (t , w ) 



d 2 G 



dw 2 



d 2 G 



dt 2 



d 2 G 



dwdt 



-i 2 



Then one checks the positivity condition [H]], in particular when the discriminant 
D (t ,w ) of the second derivatives of G at this point is positive. We can now obtain 
the approximate expression for the starting amplitude S (g,a,n) (^|): 



S Ap (g,a,n) 



where explicitly 



2vre- 1 / 2 



/ (x ,a,n) + 1 " 
D (x , a, n) 
H (x , a) exp [-/ (x Q , a, n)) 



1/2 



X 



n+1 



(20) 



D (x , a, n) 



~Xn 



df (x , a, n) 



+ f(x ,a,n) 



+f (x , a, n) [1 + / (x , a, n)) 
with x (g, a, n) , and F defined above in eq. ( |TT| ) and respectively. 



(21) 
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Notice that the expression (E^) is valid for functions H (x, a) that fulfill the following 
general conditions: 

1) the relation (|T7D must be biunivocal. 

2) D (x a , a, n) has to be positive and x a d f( x £> a ' n ) — 2 f (x a , a, n) has to be negative 
in x . 

3) / (x , a, n) ^> 1. This condition is fulfilled in particular when n ^> 1. 

These conditions provide the range of values of the involved parameters where the 
approximate solution (BO) is valid. Then, the calculation of the GBT simply consists 
on solving the implicit equation fll7j ) to obtain the saddle point expression. In general, 
the functional complexity of H can add constraints on the range of the parameters. 

We are now prepared to study the amplitude coming from the delta potential, 
namely the wave packet 

$(Q) = F ^ + 1 f°° ex P (-Qx)dx 



Q Q Jo l + vx/2[l-exp(-2/xm)} 

POO 

ee / H {x, v, m) exp (— Qx) dx (22) 
Jo 

On the other hand, the GBT is able to deal with 

/>oo 

$n (Q) = / x n H (x, v, m) exp {-Qx) dx (23) 
Jo 

The relation between (|22|) and (p3|) 



*»w) = (-r^w) 



can be inverted obtaining 



/» n—l 
dQ--- dQ1> n (Q) + J2c P (v,m)Q p (24) 

v ^J. / P=0 

n 

where the finite sum come from the indefinite integrations. Notice that all the co- 
efficients vanish whenever the Laplace transform (|22"D fulfill the following asymptotic 
behavior 

lim $ (Q) = 

Q— >oo 

Moreover, the expression ( ^4|) is valid for any value of n, in particular when n ^> 1. 
Consequently, the approximate solution reads 

$4 P (Q) ~ (-)" JdQ---fdQ*i? (Q) (25) 
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where, for n ^> 1 



$ n (Q) ~ $f (g) 



27r/e [j (x , n) + 1\ , n+ i , , , , . 

— r (x ) H (x , v, m) exp -/ (x a , n)\ (26) 

D(x ,n) 

is the approximate solution provided by GBT, with the D and / functions defined 
above flf§) and (0). 

For n sufficiently large and Q > 0, the explicit expression of the saddle point be- 
comes x = {n + 3/2) /Q. Then, we can substitute this expression into (|6|), obtaining 
the following approximate solution for the expression 



A V^(n + I)"/ 2 (2 + </ 2+1 / 2 exp (-n - 3/2) 

n W) Q n + 1 {l + v(n + 3/2)/2Q[l-exp(-2Q/(n + 3/2)m)}} 
T(n + 1) 1 
Q"+! 1 + G (n, Q) 

If < g < oo, then < G (n, Q) < v/m < 1 and we can expand to obtain 

(Q) - r (n + 1) f; [ ~ Hn ^ )/2]P £ (f) R* exp (-2*Q/ (n + 3/2) m) 
p=o ^ k=0 ^ ' 

(27) 

It is illustrative to separate the perturbative and nonperturbative contributions to 
this expression. To this end we analyze the first contribution, e.g., the term k = 0. 
This correspond to the asymptotic solution on Q . After, we add the terms with k ^ 0, 
corresponding to nonperturbative corrections. 

For k = 0, by solving the n-integrations, we obtain 

qGbt (0) „ ,^ [-v(n + m/2] P T( P + l)T(n + l) ^ 1 ^ [-,/2f 
° QP+1 r(n + P+ l) -Q^ Q Q P L[P+l) 

where in the last step, we have explicitly taken the limit n — > oo. Then, from eq.(l22f) 

oo 1 

^ G o BT (Q) * £ ("« W))" r (P + X ) = = TTyT (1/a (g)) Ei (1, l/ a (Q)) 

Consequently, we have recovered the approximate solution (|7]) provided by the 
conventional Borel transform in the region where this one provides sensible results. 
To determine the nonperturbative corrections, we must solve the n— integrals 

I n+p+l (\a\,Q) = JdQ... JdQ eXP g ( ~y } ; a^O 



to get 



In+p+i (\a\ , Q) 



\n+p / \n 



d n+p 



T(n + p + l) \a\ n ^db n +P 



exp(|a|6) Ei (n, \a\ (b + Q))\ 



b=0 



where Ei (n, x) is defined in eq. ([8]). Then, one can carry out the (n + p) derivations 
valued at b = 

Wi (M , Q) = \a\ p+1 exp (- \a\ Q) (-) n G(n+p + l,p + 2, \a\Q) (28) 
where G is the confluent hypergeometric functions of second class. 



Hence, replacing the expression fl27D into ( ^5|) and using eq. (|28|) we can write 



GET 



(Q) 



GET 



(Q) 



oo p 



+ lim QT (n + 1) V V 

p=l fc=l 

2k -Q 



exp 



(ri + 3/2)m 



G 



n + p + l,p + 2, 



-u (n + 3/2)/2] ? 
2A; 



2A- 



(ra + 3/2)m 



P+i 



n + 3/2) 



To perform the limit n — > oo, we have used here the integral representation of G . 
Then, we conclude that the approximate solution reads 



GET 



(Q) 



(Q) 



p=l k=l 



\k+p 



1 P+l 



tf p+1 I 2 




(29) 



where K p is the Bessel function of second class. 



The approximate solution ( p9|) converges rapidly to the exact solution (gzj). The 
nonperturbative corrections recover the corresponding dependence on m . The devia- 
tion with respect to the exact solution obtained is certainly small within all the range 
of Q. In particular, the second order approximation suffer a deviation lower than 4% 
f° r Q ^ 2 (see figure 1). 

The test of the accuracy provided by the GBT is presented in Figure 1 where we 
have compared the ratio Fg BT (Q) / F (Q) (line c) between our approximate analyti- 
cal solution (p9|) up to second order ( k — 2 ) and the numerical integration of the 
exact expression (p2|) , with the corresponding to the Borel resummation technique 
Fg s (Q) /F(Q) (line a) and optimization of the PT F M=4 (Q) /F(Q) (line b). This 
comparison has been performed for the particular values of the parameters m = 3 and 
v = 0.1. 

Let us now turn to the case of an attractive potential. This corresponds to v < 0, 
and can be studied by changing v — > —\v\ in the previous definition of the potential. 
In this case the expression for the wave packet (||) is not Borel summable because now 
its conventional Borel transform (ra) has a pole that does not allow the integration. 
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This situation is similar to the one presents in infrared perturbative QCD |T^[T3|[[n[ . 



In this case the renormalon technique usually extends the result of Borel summation 
(0) by using the principal value (PV) prescription. Then adds a term, coming from the 
ambiguity generated by the pole, that represents the non perturbative contribution. 
Consequently, the series (0) results in the following approximate solution 

*ns (Q) = n s v (Q) + Cexp (1/a) (30) 

where 

n s v (Q) = exp (-1/a (Q)) Ei (1, 1/a (Q)) 



and 



W { l,l/a { Q)) = P.V.I "fiE^ (31) 



The first term in (|30|) has the same perturbative asymptotic expansion as the exact 
function (0) and the constant C is obtained by using purely nonperturbative methods 
or is extracted from experiment. The optimal value obtained is C = —0.06 (see |]). 
However, the results ( |3"UD diverges strongly at small Q (see figure 2, line a). 

On the other hand, the alternative approach of the modified perturbation theory 
obtain the best approximate results for the value of the parameter fi = 3.6. The 
obtained result presents an improvement with respect to the previous approximate 
solution reducing the deviation at small Q to approximately 20% but its lose the 
asymptotic exact behavior (see figure 2, line b). 

The GBT follows the same steps presented in the case of the repulsive potential. 
Notice that previously, the term k = generates the asymptotic solution series but 
now this series in non alternate. In this regime the approximate value for this sum 
is given by expression ([H]). The different behavior is related to the non perturbative 
contribution. As it was mentioned above, the renormalon technique shows that this 
contribution can be estimated through the ambiguity generated by the non analyticity 
of the conventional Borel transform. However, as we have shown in the case of a 
repulsive potential, the origin is in the remaining terms with k ^ in the expansion 
of the binomial expression (|27|). In fact, those terms can not be taking into account in 
the renormalon technique. Thus, our approximate solution reads 



n °° p / \ 

*snQ)-nv«)-4i:i:m 

1 1 P =i fc=i \ / 



2m 
k 



Q 



p+i 



As is was the case for the repulsive potential, our approximate solution has the 
correct behavior for all the range of the momentum Q. In fact, the deviation, when 
one considers up to second order (k = 2), is lower than 3.5% for Q > 2 (see figure 2, 
line c). 
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In Figure 2 we present the corresponding comparative plots of the ratio between 
the approximate solutions Fj!f s , ^=3.6, Fns T an d the exact one F for the attractive 
potential. 

In Summary, we have presented an analytic expression for the wave packet valid for 
all the range of Q, based on the GBT. In the case of the summable Borel theory, our 
results show explicitly how the GBT complete the perturbative solution provided by 
conventional Borel transform by means of the adequate incorporation of nonperturba- 
tive contributions. In the non summable Borel theory case, besides obtaining sensible 
results, we have shown that the real origin of the nonperturvative contributions is not 
the Borel ambiguity. 
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Figure 1: numerical comparation for Borel summable solutions 



Figure 2: numerical comparation for no Borel summable solutions 
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